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Abstract 

The COBRA collaboration operates a demonstrator setup at the underground facility LNGS (Laboratori Nazionali 
del Gran Sasso, located in Italy) to prove the technological capabilities of this concept for the search for neutrinoless 
double beta-decay. The setup consists of 64 (1 x 1 x 1) cm 3 CZT detectors in CPG configuration. One purpose of 
this demonstrator is to test if reliable long-term operation of CZT-CPG detectors in such a setup is possible. The 
demonstrator has been operated under ultra low-background conditions since more than three years and collected data 
corresponding to an exposure of 218kg-days. The presented study focuses on the long-term stability of CZT detectors 
by analyzing the intrinsic, fourfold forbidden non-unique 113 Cd single beta-decay. It can be shown that CZT detectors 
can be operated stably for long periods of time and that the 113 Cd single beta-decay can be used as an internal monitor 
of the detector performance during the runtime of the experiment. 
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1. Introduction 

Neutrino oscillation experiments have shown that neu¬ 
trinos have a finite rest mass, but are unable to put a 
number on how massive they are. The COBRA (Cadmium 
zink telluride 0-neutrino double Beta-decay Research Ap¬ 
paratus, [I]) collaboration searches for neutrinoless dou¬ 
ble beta-decay (Ozx/3/3) to determine the effective Majorana 
mass of the electron neutrino by a measurement of the half- 
life of the 116 Cd decay using CZT (cadmium zinc telluride) 
semiconductor detectors operated at room temperature. 
The operation of CZT detectors in a CPG (coplanar grid) 
configuration [2] enables to use large-volume detectors that 
contain the isotope of interest and hence results in large in¬ 
trinsic detection efficiencies. As Oidecay is a very rare 
process with a half-life greater than lx 10 25 years in case of 
116 Cd, the setup has to be an ultra low-background exper¬ 
iment with a large installed detector mass and which has 
to be operated for several years to acquire the needed ex¬ 
posure. To reach the required sensitivity, it is necessary to 
provide extremely stable operation over many years of data 
taking. One of the main goals of the COBRA demonstra¬ 
tor is to test if a stable long-term operation of CZT-CPG 
detectors is possible. The stability is tested based on a 
dataset of 218kg-days exposure, acquired between Oct’ll 
and Feb’15. 


* Corresponding author 

Email address: daniel.gehre@tu-dresden.de (D. Gehre) 


2. The COBRA demonstrator 

The COBRA demonstrator [3] consists of 64 CZT-CPG 
detectors, each with a volume of 1 cm 3 arranged in four lay¬ 
ers of 4x4 detectors. The four layers are referred to as LI, 
L2, L3 and L4 in the following. The total installed detector 
mass is 380 grams. The continuous operation of the detec¬ 
tors under a controlled, ultra low-background environment 
allows for an identification of the different background con¬ 
tributions as well as to develop methods to discriminate 
such background events from the signal searched for. The 
detectors are installed in a (60 x60 x60) cm 3 combined low 
activity lead/copper castle and are kept in a nitrogen at¬ 
mosphere at room temperature. Furthermore, each detec¬ 
tor is permanently powered by its respective bias voltage 
(BV) and grid bias (GB). 

3. Measurement strategy 

Since the detectors are made from natural Cadmium 
containing 12.2% 113 Cd and as they are operated in an 
ultra low-background environment, the intrinsic, fourfold 
forbidden non-unique single beta-decay of 113 Cd( 3+ ) —> 
113 InU+) is by far the strongest signal. The decay of 
113 Cd is a direct beta-decay into the ground state of 113 In. 
Hence, the spectrum is continuously spread over the al¬ 
lowed energy range from zero to the (J-value of the decay, 
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Q = 322.2 ±1.2 keV [2j. It can be used to monitor the sta¬ 
bility of the detectors by analyzing the decay rate of this 
isotope. The theoretically expected rate of 113 Cd decays 
is about 400 per day and detector (assuming 5.9 g detec¬ 
tor mass, 5.0 at% zinc concentration and 100% detection 
efficiency). Currently, the shape of the 113 Cd spectrum 
in the range below 100 keV is not well known and experi¬ 
mental data is strongly requested [5H7]. The long half-life 
of 113 Cd of (8.00 ± 0.35) xlO 15 years [3] and its homoge¬ 
neous distribution inside the CZT detector should result in 
a constant decay rate over the time of operation. Changes 
in the measured decay rate can indicate an alteration of 
the detector properties which could have an affect on the 
overall performance of the experiment. 

4. Selection criteria and measurements 

Due to the uncertainty in the shape of the 113 Cd spec¬ 
trum, the stability analysis is limited to runs with the 
lowest common energy threshold for each detector. The 
energy threshold defines the amplitude of the pulse shape 
that must be exceeded to trigger the data acquisition sys¬ 
tem for storing events. The thresholds for the detectors 
have been adjusted on a run-by-run basis. This is done 
to reach the lowest possible energy threshold. Thus, the 
available lifetime depends on the observed energy inter¬ 
val for each detector. Runs with thresholds above 250 keV 
are discarded for the stability analysis to avoid too large 
statistical uncertainties of the detected count rate. An 
optimization procedure is applied to identify the optimal 
threshold -Eopt to maximize the count rate and to min¬ 
imize the statistical uncertainties. In this approach, the 
total number of counts are maximized for each detector as 
a function of the available lifetime per energy bin. 

In Figure [l] the total spectrum (blue) of one detec¬ 
tor is shown. The available, threshold-dependent detector 
lifetime in days is plotted in red. Data taken below the 
optimal threshold are not used for the stability analysis. 
The lifetime distribution shows that for lower energies the 
available lifetime of the detector drops fast and for ener¬ 
gies below 60keV it is less than 20 days. For the stability 
analysis, an additional ±20keV offset is applied to A 0 pt- 
This accounts for the limited energy resolution in that en¬ 
ergy range. The remaining part of the spectrum (green 
histogram) is used for the stability analysis. Runs with 
energy thresholds higher than A op t are discarded. 

A further limitation is caused by the working principle of 
the coplanar grid. The pulse-shape analysis used in the 
COBRA demonstrator allows for the reconstruction of the 
depth of interaction (DOI) as well as to judge if the in¬ 
teraction happened at the lateral surfaces of the detector 
BM- It is known that in the vicinity of the anodes the 
presence of the grid bias field distorts the movement of 
the charge carriers. This distortion alters the energy- and 
depth-reconstruction for energy deposits close to the an¬ 
odes. For this reason, a second cut (static depth-cut) is 
used to limit the fiducial volume to the section where the 
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Figure 1: Exemplarily, in blue the total 113 Cd spectrum of one detec¬ 
tor is plotted. The green part of the spectrum is used for the stability 
analysis. The threshold dependent lifetime of the detector is shown 
in red. The different lifetimes per energy bins are caused by the 
operation of the detector at adjusted thresholds in the subsequent 
runs. 


energy and depth are reconstructed correctly. The static 
depth-cut selects the region >0.2 and <0.97 of normalized 
depth and, hence, excludes the effect of the hole shift [8j 
in the vicinity of the anodes. 

5. Observed 113 Cd decay rate variations 

The operation under ultra low-background conditions 
and the analysis of the high-energy part of the spectrum 
allows to neglect further contributions from other back¬ 
ground sources in the analyzed energy range from 100 keV 
to 350 keV for the stability analysis. It is therefore as¬ 
sumed, that the 113 Cd signal region is not contaminated 
by background. To analyze for the stability of the detected 
113 Cd decay rate, the data is partitioned into ten-days real¬ 
time intervals. The lifetime of the detectors has to be at 
least five days for the respective interval to be considered. 
For each detector, the optimal threshold-cut and the static 
depth-cut are applied. The remaining 113 Cd counts are 
summed up and normalized to the lifetime in the respec¬ 
tive time period. Due to the different optimal thresholds 
of each detector, the observed absolute decay-rate deviates 
from one device to another. 

The measured 113 Cd decay rates of detectors Det4, 
Det5, Det6 and Det8 of LI are shown in Figure [2] (red 
marker). The detectors were operational at all times. Miss¬ 
ing data points indicate an operation above the optimal 
energy threshold in the specific time period which led to 
a removal of the run. A robust, linear approximation is 
superimposed. The grey boxes indicate excluded time pe¬ 
riods (Jul’13 - Nov’13 and Jul’14 - Sep’14), which were 
affected by electromagnetic interferences that caused an 
increased baseline noise. This noise led to a higher re¬ 
jection rate of triggered events and, thus, accidentally re¬ 
moved low-energetic physics events. The 113 Cd rate is es- 
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Figure 2: The threshold dependent 113 Cd decay rate of Det4, Det5, 
Det6 and Det8 of LI. The data is partitioned into ten day intervals, 
whereas the rate is calculated in counts/day in the analyzed part of 
the spectrum and varies due to the different optimal thresholds. The 
uncertainties on the rate are purely statistical in nature. A linear 
approximation is superimposed in blue. The grey boxes indicate 
excluded time periods. 


pecially affected since the deposited energy is always below 
322 keV. Hence, it was decided to exclude the data from 
these time periods for the evaluation. 

Det4 shows the strongest decrease of the 113 Cd rate in 
LI (Aj ) et 4 = -5.2% per year), whereas Det6 is the only 
detector that shows a slight increase of the count rate 
(A Det6 = + 1.4% per year). All other detectors are in be¬ 
tween those values for the observed rate variations. For 
Det8, the measured decay rate of the 113 Cd is basically 
constant over the full time period (excluding the electron¬ 
ically disturbed runs). This detector yields one of the 
largest data sets with more than 810 days accepted life¬ 
time. As the detectors are produced using material from 
different batches, the single layers have to be considered in¬ 
dependently. The layers were installed subsequently with 
several months (LI and L2) or more than one year delay 
(L3 and L4). Furthermore, L3 is still not performing as 
well as the other detector layers and has been excluded 
from the following considerations. This layer is generally 
more affected by electronic noise and the overall detector 
performance lags behind those of the other detectors in 
terms of energy resolution. Additionally, the considered 
exposure of L3 and L4 is only about one third of those of 
LI and L2, which results in fewer time periods and larger 
uncertainties for the linear approximation. For LI, the av¬ 
erage accepted lifetime of the detectors is roughly 700 days 
whereas it is for L2 630 days, for L3 300 days and for L4 
it is only 270 days. Nevertheless, the detected count rate 
variations are relatively low. Beside the excluded time pe¬ 
riods no significant outliers are found that would indicate 
a drastic change of the detector properties or of the data 
acquisition chain. This can be seen in Table [lj 

The relative change in count rate per year and detector 
is plotted in Figure [3] The error bars indicate the uncer¬ 
tainty of the linear approximation on the data which is 


Table 1: Relative 113 Cd decay rate variations per year for each de¬ 
tector. LI, L2 and L4 are considered for the stability analysis. L3 
is excluded due to general performance problems in the low-energy 
range. 


LI 

Rel.rate change/yr 

L2 

Rel.rate change/yr 

Detl 

1.000T0.004 

Detl7 

1 . 001 ± 0.001 

Det2 

0.977±0.016 

Detl8 

l.OlOiO.OlO 

Det3 

0.960T0.019 

Detl9 

0.993T0.001 

Det4 

0.948T0.052 

Det20 

1.006T0.006 

Det5 

0.976T0.024 

Det21 

1.025±0.025 

Det6 

1.014T0.015 

Det22 

0.973T0.016 

Det7 

0.996T0.004 

Det23 

0.967T0.018 

Det8 

1.000T0.002 

Det24 

1.015T0.015 

Det9 

1.005±0.005 

Det25 

0.991T0.012 

DetlO 

0.998T0.005 

Det26 

1.035T0.032 

Detll 

0.984T0.016 

Det27 

1.003T0.007 

Detl 2 

1.004T0.004 

Det28 

1.033T0.033 

Detl 3 

0.999T0.001 

Det29 

1.031T0.029 

Detl4 

1.001T0.003 

Det30 

0.969T0.031 

Detl 5 

0.998T0.002 

Det31 

0.961T0.022 

Detl 6 

0.943T0.057 

Det32 

1.056T0.056 

L3 

Rel.rate change/yr 

L4 

Rel.rate change/yr 

Det33 

0.982T0.001 

Det49 

0.980T0.020 

Det34 

1.004T0.014 

Det50 

0.997±0.012 

Det35 

0.966T0.034 

Det51 

0.980T0.020 

Det36 

0.695T0.305 

Det52 

1.030±0.030 

Det37 

0.976T0.024 

Det53 

1.032±0.032 

Det38 

0.891T0.109 

Det54 

0.985T0.006 

Det39 

0.892T0.108 

Det55 

0.951T0.006 

Det40 

0.928T0.072 

Det56 

1.042T0.031 

Det41 

0.960T0.010 

Det57 

0.951T0.049 

Det42 

0.849T0.070 

Det58 

0.983T0.002 

Det43 

0.916T0.035 

Det59 

0.908T0.032 

Det44 

0.964T0.004 

Det60 

1.022T0.022 

Det45 

0.839T0.161 

Det61 

0.978T0.002 

Det46 

0.857T0.071 

Det62 

0.989T0.011 

Det47 

0.705T0.295 

Det63 

1.012T0.012 

Det48 

0.733T0.267 

Det64 

1.038T0.038 


basically related to statistical uncertainties. 45 out of 48 
detectors from LI, L2 and L4 show a change in the decay 
rate well below ±6% per year. The shorter overall lifetime 
and the higher optimal thresholds of the detectors of L4 
affects the accuracy of the approximation as only a smaller 
part of the spectrum is analyzed and fewer accepted run 
periods are available. Nevertheless, a weighted mean of 
pi = (0.995±0.004) can be reported, which indicates an 
almost perfect stable operation of the detectors. 

Figure [4] displays the lifetime-weighted relative change 
rate distribution for all layers (blue) and for LI, L2 and 
L4 in red. The distribution is clustering around one with 
a slight shift towards values less than one. 

To analyze for the cause of the observed changes, the 
spectral shape and the depth distribution of different time 
periods of the experiment are compared. To accomplish 
this, the number of available intervals is sectioned into 
three parts: the start, middle and end of the lifetime. For 
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Figure 3: The relative 113 Cd count rate change on a per year basis is 
plotted in blue. In red the respective optimal threshold per detector 
is shown. 



Figure 4: The lifetime-weighted distribution of the relative rate 
change of all Layers (blue) and Layerl, Layer2 and Layer4 in red. 


the three time periods, the counts per keV are summed 
up bin-wise and are normalized to the respective lifetime 
per bin for each detector. In Figure [5j the three spectra of 
the detectors are plotted in red (first third), green (second 
third) and blue (last third). In grey, the part of the spec¬ 
trum that lies above the (F O pt+20) keV is marked. Only 
this part is used for the stability analysis. In the presented 
cases the red spectra (first third of the data taking period) 
tends to he slightly above the green and blue ones - but 
no clear tendency is observable if one takes the other de¬ 
tectors into account. In general, no shift or deformation of 
the spectral shape has been found comparing the different 
time periods. 

For the stability analysis of the reconstructed depth 
distribution, the same time sectioning and least common 
energy threshold-cuts are applied. In contrast to the decay- 
rate based stability analysis that implies the static depth- 
cut, no depth dependent cuts are applied here. This al- 



Figure 5: Comparison of the spectral shapes of the 113 Cd decay 
for different time sections. The data is partitioned into three time 
periods: first third (red), second third (green) and last third (blue). 
The grey part of the spectrum has been used for the stability analysis. 
The lower inset shows the relative deviation of the middle (green) 
and the end period (blue) with respect to the start period (red). 


lows for an analysis of the full depth distribution of the 
recorded interactions. No time-correlated impact on the 
reconstructed depth is visible as can be seen in Figure [6j 
The fact that the shape of the depth distribution looks dif¬ 
ferent for the single detectors is caused by their different 
thresholds. The lower the threshold, the higher the total 
count rate. The asymmetry of the reconstructed depth, 
which shows a shift of the reconstructed events towards 
the cathode, is caused by imprecisions of the reconstruc¬ 
tion algorithm for lower energies. 

6. Conclusion 

The presented study is based on the analysis of the in¬ 
trinsic, fourfold forbidden non-unique single beta-decay of 
113 Cd with a Q-value of (J=(322.2 ± 1.2) keV. The result 
obtained is limited to the energy range from approximately 
100 keV to 350 keV. The analyzed dataset represents a life¬ 
time of roughly 3.5 years with a total acquired exposure of 
218kg-days, which is the longest measurement ever taken 
with CZT detectors under these conditions. 

Based on this dataset, count rate variations in the lower 
percent-per-year range ( A mean = (—0.5 ± 0.4) % per year) 
are found for 45 out of 48 detectors. The majority of the 
detectors are working without any glitches or sudden per¬ 
formance changes. The reproduction of the spectral infor¬ 
mation as well as the depth reconstruction of the detectors 
is not affected and is stable over the analyzed time period. 
The 16 detectors of layer L3 are excluded due to known 
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Figure 6: Comparison of the depth distribution of the 113 Cd decay 
for different time sections (anodes at ’O’, cathode at ’1’). For the 
stability analysis the static depth cut has been applied. 


performance issues. Disregarding L3, the overall perfor¬ 
mance and stability of the setup is excellent. Currently 
it is not clear what causes the observed minor alterations 
of the count rate and if this is seen for higher energies as 
well. To analyze for efficiency alterations in higher energy 
ranges, a continuous monitoring of the detectors based on 
precisely positioned calibration sources is required. Nev¬ 
ertheless, under the precondition of ultra low-background 
operation, the analysis of the intrinsic 113 Cd decay offers 
a unique way to monitor the detector performance during 
the runtime of the experiment. 

Finally, it can be concluded that CZT detectors can be op¬ 
erated stably under ultra low-background conditions over 
time scales of at least several years. This result shows that 
it is possible to move towards a large-scale experiment as 
proposed by the COBRA collaboration. 
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